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THE FITNESS OF THE STRONGEST INDIVIDUAL IN 
THE SUBCRITICAL GMS MODEL 


CAROLINA GREJO, FABIO MACHADO, 

AND ALEJANDRO ROLDAN-CORREA 

Abstract. We derive the strongest individual fitness distribu¬ 
tion on a variation for a species survival model proposed by Guiol, 
Machado and Schinazi We point out to the fact that this dis¬ 
tribution relies on the Gauss hypergeometric function and when 
p = I on the hypergeometric function type I distribution. 


1. Introduction 

We consider a discrete time model beginning from an empty set. At 
each time n > 1, a new species is born with probability p or there is 
a death (if the system is not empty) with probability q = 1 — p. Let 
Xn be the total nnmber of species at time n. Xn is a random walk 
on that jnmps to right with probability p and jnmps to left with 
probability q. When Xn is at 0 the process jnmps to 1 with probability 
p or stays at 0 with probability 1 — p. We assign a random nnmber 
to each new species. This nnmber has a nniform distribntion on [0,1]. 
We think of this nnmber as a htness associated to each species. These 
random nnmbers are independent to each other. When a death occnrs, 
the individnal with lowest htness dies. This model, latter denominated 
GMS model, was hrst proposed and stndied in Gniol ef a/ [5]. Some 
interesting variations were fnrther stndied in Gniol ef a/ [B], Ben Ari et 
al |2] and Skevi and Volkov [TO] . 
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In Guiol et al [5j it is shown that there is a sharp phase transition for 
p > 1/2. For Rn, the set of species with htness higher than fc = at 
time n approachs an uniform distribution in the following sense. For 

fc<a<b<l 

\Rr,n{a,b)\ 

n^oo Tl 

On the other hand every specie born with htness less than disappear 
after a hnite (random) time. The set of species present in the system 
whose htness is smaller than fc becomes empty inhnitely many times. 

Here we focus on the case p < 1/2 in order to understand better 
the dynamics of this model. In this case, the process is recurrent 
and the system becomes empty inhnitely many times. Therefore it is 
not interesting to study the distribution of the htness of the species 
which are alive on the system in the long run. An interesting point is 
to study the distribution of the htness of the strongest individual on 
each excursion between the epochs when the system becomes empty. 

We propose a variation for the GMS model by considering that each 
time the system becomes empty, a set of m individuals are introduced 
with independent set of htness. This variation is meant to reinforce 
competition among species before the system becomes empty again. 

2. Results 

We deduce explicitly the distribution of the htness of the strongest 
individual on excursions between the epochs when the system becomes 
empty. The last individual to die before the system becomes empty is 
the strongest on that excursion because the hrst ones to die are those 
individuals with the smallest htness. 

Observe that some excursions may have length 2. When this hap¬ 
pens, the individual who is born, dies right away without competing 
with any other individual. To ensure that each excursion has com¬ 
petition among individuals in a sort of natural selection process, we 
introduce a change-over on the model; Each time after the system be¬ 
comes empty, m independent new species are placed on the system 


p{b — a) a.s. 
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(a) GMS(1) 


(b) GMS(IO) 


Figure 1. Histograms of the fitnesses of the strongest 
individual on GMS(m) after 200,000 births and deaths 
for p = 1/4. 


(instead of just 1) with probability p, or the system stays empty with 
probability 1—p. We denote this variation by GMS(m). In this set up 
GMS(l) is the original model. 


Figures Im and |lb| show the role of the competition on the distribu¬ 
tion of the fitness of the strongest individual on each excursion. Short 
excursions are more commom on GMS(l) than on GMS(IO). That be¬ 
haviour favors individuals with lower fitnesses to be the strongest ones. 
Gompetition introduced in GMS(IO) avoids that. 

The next result computes the htness distribution of the strongest 
individual to die right before the system becomes empty on GMS(m) 
model. It is shown in terms of the hypergeometric function of Gauss 
(see Luke 0). This function is denoted by 2 Fi{a, b; c; z), namely, 

ia)k{b)k 


2^1 (a, 6; c; z) = ^ 


k>0 


(c)fc k\ ’ 




( 2 . 1 ) 


where a, b, c, are real numbers with c ^ 0, —1, —2,..., and (a)k is the 
coefficient Pochhammer, namely. 


(a)k = a{a + 1) ■ ■ ■ (a + k - 1) (a)o = 1. 

Theorem 2.1. Let p < 1/2 and Zm be the fitness of the strongest 
individual before the system becomes empty on GMS(m) model. Then 
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Zra is a random variable with distribution 

( ifT) ‘YY! _j_ 1 \ 

—, —;m + l]Apqt\ 


0 < t < 1. 


Corollary 2.2. Let p < 1/2 and Z be the fitness of the strongest 
individual before the system becomes empty on GMS(1) model. Then 

1 — ^/T^^Apqi 


P[Z <t] = 


2p 


0<t <1. 


For p = 1/2, Z follows a Beta distribution i?(l,l/2). 

Remark 2.3. By Tlieorem|2.1|we have Zm density probability function 


IS 


dt 


{qtT 2 Fi ( m + 1; Apqt 


rr, rv,_1 ^ f VO m + 1 

= mq t zFi I —, ^ ;m;4pqt 


( 2 . 2 ) 


where the last line have been obtained by using Abramowitz and Stegun 
[U Eq. 15.2.4], When p = g = 1/2, the distribution of 1 — Z^ is known 
as hypergeometric function type I distribution (see Gupta and Nagar 
m p. 298]). 

( TTl 7Tl 1 

Corollary 2.4. E[Zm] = 1 - — 2 F 1 ( —, —-—;m + l;4pg 

m +1 \ 2 2 


Remark 2.5. Considering Corollary 2.4 when p = 1/2, by using the 
following equality (see Gradshteyn and Ryzhik |U Eq. 7.512.11]) 

r(7)r(7-«-/d) 


2^1(0,/ 3 ; 7; 1) = 


r(7-«)r(7-/^) 


we have that 


^ fm m + 1 

E Zm = 1-2^11 —,-: m + 2; 1 

^ ^ m +1 V 2 ’ 2 ’ ’ 


= 1 - 


2^{m + 1) 


r(m + 2)r(3/2) 


r(^ + 2)r(^^) 


= 1 


^/^T m\ 


2"^+ir(^ + l)r(^) 
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where the last line has been obtained by using r(3/2) = a/t/2. Now, 
using the duplication formula, namely. 


we get 


T(2z) 


m + Drw 

2l-2zy^ 


E[Z^] = 1 


2 

{m + l)(m + 2) 


3. Proofs 


Proof Theorem 2.1. For n 


0,1,... we dehne 


Tn = inf{/ > 1 : Xn+i = 0, = 0}. 


In words is the length of a excursion from 0 to 0. As the process 
Xn is homogeneous, the distribution of does not depend on n so we 
consider the random variable r := tq. Besides, as p < 1/2 we have that 
P[r < cx)] = 1 and 

p[r = k+1] = F[T_m = ^] = J k>m, k+m even, 

where T_m is the time of the first visit to —m for a random walk on Z 
beginning at 0. (See Bhattacharya and Waymire |3]) 


If r = A; + 1, we see along that excursion, extra births and 
deaths. The last death corresponds to the individual with the strongest 
htness among all that were born. Hence, 


P[-^m < A] = P[r = k + l]P[max(Yi,..., Ym+k) < t], 

k=m 

where Ti,.... Ym+k are i.i.d. uniform random variables on [0,1]. There¬ 
fore, 
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nZm <t] = (kL even} 


k=m 


E UL 

9 / - 


m [21 —m\ It,, 

‘ 'p qh^ (k + m = 2[l>m) 


l=m 


21 — m\ I 


m /m + 2A , 

m + 2] \ m + j J 


{qtr±^^{^^^^\{pqty 
^ m + 2j V + J 


{qtT 


{jn)2k {pqty 


= (gtr2Fi(f,^;m + l;4pgt) 


where the last line has been obtained by using ( 0 ) 2 ^ = (|)^ 

and the definition of Gauss hypergeometric function. □ 


Proof Corollary \2.S\ It is a particular case of Theorem 2T when m = 1. 
In this situation 


qt 2 F 1 1 -,l;2;4pgf 


qtY 


(l/2)fc(l)fc {ipqty 


k>0 


(2). 

{2k)\ 


-T 

Pk^o (^ + 1)' 


k\ 

{pqt)’^+^ 


^/l - Apqt 

2p 


where the last line has been obtained by using (l)fc(l/2)fc = 2 ^^(2fc)! 
and the result given in Prudnikov et al [9l Eq. 5.2.13.8]. □ 
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Proof Corollary S.fl 


¥\Zm] = Jq P[Zm > t] dt 


m m+l . 
2 ’ 2 ’ 


'^^;m + l;4:pqt) dt 
m + 2] 4pg) 



m m+l . 
2 ’ 2 ’ 


where the last line has been obtained by nsing the resnlt given in Grad- 


shteyn and Ryzhik |U Eq. 7.512.11], 


□ 


Acknowledgements: The anthors are thankful to Daniel Valesin and 
Rinaldo Schinazi for helpful discussions about the model. Thanks are 
also due to the anonymous referee for his/her constructive comments, 
leading to an improved presentation. 


References 


[1] M.Abramowitz and I.Stegun.Handbook of Mathematical Functions. National 
Bureau of Standards. (1970). 

[2] I.Ben-Ari, A.Matzavinos and A.Roitershtein. On a species survival model. Elec¬ 
tronic Communications in Probability, 16 , 226-233. (2011). 

[3] R.N.Bhattacharya and E.C.Waymire. Stochastic Processes With Applications. 
Society for Industrial & Applied, New York. (2009). 

[4] I.S.Gradshteyn and I.M.Ryzhik. Table of integrals, series, and products, sev¬ 
enth edition. Elsevier/Academic Press, Amsterdam. Translated from the Rus¬ 
sian, Translation edited and with a preface by Alan Jeffrey and Daniel Zwill- 
inger. (2007). 

[5] H.Guiol, F.Machado and R.Schinazi. A Stochastic model of evolution. Markou 
Process. Related Fields, 17 , no. 2, 253-258. (2011). 

[6] H.Guiol, F.Machado and R.Schinazi. On a link between a species survival 
time in an evolution model and the Bessel distributions. Brazilian Journal of 
Probability and Statistics, 27 , no. 2, 201-209. (2013). 

[7] A.K.Gupta and D.K.Nagar. Matrix Variate Distributions. Chapman & 
Hall/CRC, Boca Raton (2000). 

[8] Y.L.Luke.The Special Functions and Their Approximations.Academic Press, 
New Yorkvol. 1. (1969) 

[9] A.P.Prudnikov, Yu.A.Brychkov and O.I.Marichev Integrals and Series. Vol¬ 
ume 1. Elementary Functions. Taylor & Francis, London. Translated from the 
Russian by N.M. Queen (2002). 

[10] M.Skevi and S.E.Volkov. On the generalization of the GMS evolutionary model. 
Markov Process, and Related Fields, 18 , no. 2, 311-322. (2012). 








THE FITNESS OF THE STRONGEST INDIVIDUAL 


(F. Machado, A. Roldan, C. Grejo) Statistics Department, Institute of 
Mathematics and Statistics, University of Sao Paulo, CEP 05508-090, 
Sao Paulo, SP, Brazil. 

E-mail address: carolina@ime.usp.br 

E-mail address: fmachado@ime.usp.br 

E-mail address: roldan.alejo@gmail.com 



